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Abstract: Contemporary experiments in cavity quantum electrodynamics 
(cavity QED) with gas-phase neutral atoms rely increasingly on laser 
cooling and optical, magneto-optical or magnetostatic trapping methods 
to provide atomic localization with sub-micron uncertainty. Difficult 
to achieve in free space, this goal is further frustrated by atom-surface 
interactions if the desired atomic placement approaches within several 
hundred nanometers of a solid surface, as can be the case in setups 
incorporating monolithic dielectric optical resonators such as microspheres, 
microtoroids, microdisks or photonic crystal defect cavities. Typically in 
such scenarios, the smallest atom-surface separation at which the van der 
Waals interaction can be neglected is taken to be the optimal localization 
point for associated trapping schemes, but this sort of conservative strategy 
generally compromises the achievable cavity QED coupling strength. Here 
we suggest a new approach to the design of optical dipole traps for atom 
confinement near surfaces that exploits strong surface interactions, rather 
than avoiding them, and present the results of a numerical study based on 
■'^K atoms and indium tin oxide (ITO). Our theoretical framework points 
to the possibility of utilizing nanopatterning methods to engineer novel 
modifications of atom-surface interactions. 
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1. Introduction 

Over the past decade there have been numerous experimental and theoretical studies in atomic, 
molecular and optical physics HI |2] [3] |4] |5] 16] and in quantum information science Q [8] |9l 
focusing on the interaction of gas-phase neutral atoms with static or electromagnetic fields pro- 
jected by planar micro- or nanofabricated structures. The tight dimensional control that can be 
achieved via modern lithographic techniques generally enables new regimes of performance 
and scalability in integrated atomic devices, which hopefully can be used to further the success 
of recent technology breakthroughs such as chip-scale atomic clocks |ll!l and magnetome- 
ters ifTOll . 

In many such scenarios — in particular those exploiting strong coupling of trapped atoms 
with optical fields confined in photonic waveguides or resonators — robust and accurate local- 
ization of atoms near solid surfaces emerges as a critical requirement. Indeed, the prospect of 
tight control over atomic ensembles alone has driven much of the development of integrated 
atomic platforms. Often atomic localization is achieved by "scaling down" trapping techniques 
first developed in free-space, as in the manipulation of neutral atoms by magnetic field gra- 
dients created by current-carrying microwires Q or by optical dipole forces induced by the 
evanescent tails of confined laser fields lfT2l . Robust and conceptually simple, these techniques 
typically improve in strength and efficiency as the atoms are brought closer to the defining struc- 
ture(s). However, bulk surface effects become dominant at short distances. For small values of 
the atom-surface separation r, energy potentials for optical dipole traps and magnetostatic mi- 
crowire traps are roughly proportional to 1 — r/ro and ' , respectively, whereas van der Waals 
potentials have a much stronger variation ~ r^^ . Hence many integrated atomic systems are ex- 
pected to exist in an awkward regime where atoms are kept close to trap-defining structures to 
achieve tight potentials, but not so close that van der Waals effects become problematic. In the 
context of cavity quantum electrodynamics (cavity QED) with monolithic dielectric resonators 
this can be a particularly significant compromise, as it generally precludes realization of the 
strongest possible atom-photon coupling parameters, which would require stable placement of 
atoms nearly in contact with the resonator surface. For this reason and others IfTSl . van der 
Waals interactions can thus become a limiting factor in the achievable performance of a wide 
range of integrated atomic devices, and we are motivated to investigate new approaches to atom 
manipulation that can compete with the strong onset of common r^^ surface effects. It may be 
anticipated that such approaches will rely on the physics particular to atoms in the near field of 
bulk structures, rather than the reapplication of techniques developed in vacuo. 

A reasonable first step in this direction could be to explore the design of trapping schemes 
that attempt to exploit serendipitous aspects of the atom-surface interaction. For example, the 
demonstration of a positive van der Waals energy shift of a particular excited state of Cs over 
an oriented sapphire crystal that is orders of magnitude larger than typical, attractive van der 
Waals potentials lfT4l [Tsl [161 suggests a possible means of surface-specific manipulation of 
atoms. If the van der Waals, repulsive character of an excited atomic state may be accessed 



appropriately in the near-field of a surface, it would provide a valuable tool for the efficient 
localization of atoms very close to structures of interest. The anomalous van der Waals pertur- 
bation demonstrated in HJI arises from the off-resonant coupling of an atomic transition to a 
surface polariton resonance in the sapphire crystal, suggesting that analogous and potentially 
improved phenomena could be accessed in alternative systems chosen or even created using 
modern methods and insights from materials science and nanofabrication. 

In this article we describe and analyze (through numerical simulations) a scheme for using 
off-resonant coupling of gas-phase atoms to surface polariton resonances in a surface near field 
to produced van der Waals-enlianced optical potentials. In section |2] we summarize previous 
work that connects atom-surface van der Waals effects to a perturbative interaction with surface 
field modes. Then in section[3]we describe a general approach to harnessing these effects with 
laser fields to produce surface enhanced optical potentials. Section H] illustrates this approach 
with simulations of specific systems based on ^^K and a planar indium tin oxide (ITO) surface. 
And in concluding, we mention some possible paths toward engineering general, useful van der 
Waals interactions. 



2. Atom-surface interactions 

QED interactions between atoms and/or bulk materials have been analyzed in many contexts 
over almost 80 years (e.g. ifTTlfTSl ). but this section will largely summarize relevant results in 
work by WyUe and Sipe lfT9l l20l 1211 (whose perspective we adopt throughout the article) as 
they naturally highlight the importance of surface excitations in producing anomalous van der 
Waals effects. 

An atom near a macroscopic structure experiences effects often dominated by the perturbing 
interaction of the atomic dipole and radiation field 

Hi = -JID, (1) 

where p. is the atomic dipole operator and D is the transverse displacement field operator at the 
atom's position. Thus, just as atomic level shifts and the spontaneous decay of excited states 
arise from interactions with radiation modes in a vacuum, van der Waals effects are produced by 
similar physical processes, with contributions from a near field, planar surface. For example, 
in a purely semi-classical description, an excited atomic dipole near an interface will have a 
reaction to the reflected portion of its own induced radiation field. As the amplitude of the 
reflected field at the atom's position increases, so do these reactive, surface-induced effects. 

Without having to specify any geometry or optical properties of the system yet, the energy 
shift of atomic level a may be approximated by second order perturbation theory 

5Ea = -h-'V y p{B)- " f , (2) 

where P indicates the principal part, the summation is over field states B and and atomic 
states n, hcOj is the energy eigenvalue for un-perturbed field or atomic state /, p{B) denotes an 
assumed thermal distribution of the field, and the summation over vector components a and j3 
is implied. Also, the transition rate from atomic state a to n is given by Fermi's golden rule 
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We introduce the correlation functions 

Gapit) = iT2-\[Dait),DpiO)])e{t) 

<pit) - in-'{a\[fla{t),fipm\a)@{t), (4) 



where Da (f ) and Ha (f ) are interaction picture operators with respect to //, and 0(f) is the Heav- 
iside step function, and also specify their Fourier transforms: Gafiico) and a^p(a)). Eqs. (HJi 
describe the linear response of operator expectation {X{t)), given a perturbation proportional 
to X{0). And, from linear response theory, G^jg {(o) and a'^p {(o) are the generalized suscepti- 
biUties for the field and the atom. In particular, Ga,^{(o) may be identified with the expected 
displacement field induced by an oscillating classical dipole, and a^^ {(£)) the polarizability of 

the fl* atomic state by the field. The fluctuation-disspation theorem, a zero temperature approx- 
imation (appropriate for optical frequencies at room temperature), and the analytic properties 
of the susceptibilities allow Eqs. (|2|i and (O to be written 

5E„ = 5E:/ + dE', 

-^Aia"M^"ReG„p(a)„„)©(»„„) (5) 



= j^^l^'^lpmGap{COa„)&{COan) (6) 

where co„„ = (Oa — to,,. 

One of the advantages of this approach is that the contribution to Gap {(o) from an interface 
may be simply added that of free space by the superposition principle, and therefore the surface 
effects can be considered separately from the (formally divergent) free space effects. From here 
on, we make this separation and only consider the effects from field modes emanating from the 
surface. The first term on the RHS of Eq. (|5]l, dE^,/, may be associated with the polarization 
energy of the atom due to quantum fluctuations in these surface modes, in analogy to the Lamb 
shift in free space. As the term integrates the susceptibilities over imaginary frequencies, this 
shift is always real, usually negative and generally insensitive to resonances in the environ- 
ment, giving rise to the familiar, attractive van der Waals potential for ground state atoms. This 
is in contrast to the last term in Eq. (|5]l, 5E[,, which arises from the residue of the principal part 
integration in Eq. (|2]i. This term has none of the analytic constraints of the first term: it may be 
either positive or negative, and is highly sensitive to environmental resonances near C0a„- Due to 
the zero temperature field assumption, this residue exists only for excited atomic states, which 
may be degenerate with other states of the atom plus field. As it originates from global energy 
degeneracies, this term admits a semi-classical physical interpretation of the energy associated 
with an atomic dipole interacting with the in phase portion of its own, expected reflected field. 
Similarly, Eq. ^ may be viewed in analogy to the dephasing rate of an excited dipole from the 
interaction with the out of phase portion of the expected reflected field. 

The atomic polarizability can be shown to be 



a, 



wj3 



(^) = ^I r.2 _1 ' (7) 



while the field susceptibility may be constructed from a plane wave decomposition of the re- 
flected electromagnetic field produced by an oscillating dipole above a particular surface. To 
specify now some general aspects of the system, the interface is assumed to be translationally 
invariant in the [xy) plane (see figure [1) which gives 

(«) - ^ / ^ {SaSpR' + ;30+;3«-/?'') exp {2iW,z) (8) 



'Each n-level's contribution to SeI/ carries the sign of (B„n, as is suggested from Eq. (TJ- Thus for downward 
isitions, SE'J will be positive, but the net s 
that is moderately reflective at |£0„o| (see Eg. [18 



transitions, SEI/ will be positive, but the net shift will often be negative because SE'^ ~ —2\8E]/\ for any interface 
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Fig. 1. The assumed geometry of an oriented atomic dipole in vacuum at height z above a 
planar, transfer region interface that separates an underlying bulk substrate with dielectric 
constant ei from the vacuum. Adapted from fjll- 

where &> = oc^', K and Wo are the components of the vacuum wave vector at the transfer 
region-vacuum interface perpendicular and parallel to z, respectively, s and are the s- and 
p-polarized electric field vectors in the vacuum for upward (+) and downward (— ) propagating 
fields, are the Fresnel coefficients of the transfer region for the respective polarizations, and 
z is the height of the atom above the surface of the transfer region. 

Although surface interactions are typically considered in the near field limit, which leads to a 
dramatic simplification of Gap{(o), Eq. dS]) is especially useful for engineering considerations, 
as it directly connects the effects on the atom by the surface with an intuitive picture of the 
dipole interacting with all possible, reflected surface modes at frequency co. In particular, if the 
Fresnel coefficients are large over a large range of K at an atomic transition frequency ©„„, the 
corresponding effects described by Eq. ^ and the last term of Eq. ^ will be likewise large 
and carry the sign of the integrated Fresnel coefficients. 

The common case of an atom above a bulk substrate (no transfer region) with separations 
z ^ 0)mi may be considered by taking the near-field limit (c °o) in Eq. (O. In this limit, the 
Fresnel coefficients become JC-independent, and the integral may be simply evaluated: 



G,-{(o) = 2Gv.v(a)) = 2G„(to) 
1 £1(0)) -1 



(9) 



4z3 ei {(o) + 1 

where the final "image factor" is the familiar, "quasi-static," p-polarized Fresnel coefficient of 
the interface (and all other G^p (to) components are 0). In this limit, all surface modes coupling 
to the atom are evanescent (ReWo = 0), and therefore bound to the surface. Resonant surface 
excitations in all K modes are excited by the atom as ei (a)a„) — > — 1, leading to a divergent re- 
sponse in dE^j and Rmi- Significantly, an anomalously large, positive, van der Waals-like energy 
shift of excited atomic states may be realized by off-resonant coupling to these surface exci- 
tations in the near field as e(ft)o„) + 1 ^ 0+. An intuitive explanation for this effect is that the 
atomic dipole may interact with a near field image dipole (without retardation) that because of 
material resonances is both larger than and ~ 180°-shifted from the image dipole in an ordinary 
reflector. 

For example, consider a Drude metal with the dielectric function 

eoico) = 1 - . (10) 

where cOp is the metal's plasma frequency and 7 the electron collision rate. Using £d, the real 
part of the image factor in Eq. ^ reaches a minimum of -2^(-^ + (^)^) 



0) — ( ^ + ) ' = (Omin ■ Compared to 1 , image factor for a perfect conductor, this minimum 
coefficient has the opposite sign and is far larger in magnitude when cOp ^ 7. The correspond- 
ing imaginary part of the minimum image factor is —IminO- + V^-^)^''^, which notably also 
diverges for a "dissipation free" metal, 7^0. An excited atom with a downward transition 
frequency near (»,„,„ in the near field of this Drude metal will experience a positive energy shift 
proportional to that will dominate 5Ea-^ when fflp ^ 7. Also, the spontaneous emission rate 
of that particular transition will likewise be enhanced by a factor of order tT ' Sis^p proportional 
toz"^ 

3. Van der Waals atom mirrors 

The flexibility of 5E'^^ suggests that a resonant material interface may be used to engineer the 
energy potential of an atom in its near field. However, to do so we are required to find a suitable 
means of accessing these excited state effects, while keeping the atomic heating from sponta- 
neous emission manageable. 

Previously used to measure the attractive van der Waals potential for alkali atoms near trans- 
parent dielectrics 1221 . optical atom mirrors could provide a simple platform to study and begin 
to exploit anomalous interactions of atoms and surface polaritons. In a traditional optical atom 
mirror, an evanescent, coherent field is produced by total internal reflection of a laser beam 
from a planar, vacuum-dielectric interface (see figure |2]l. If the frequency of the laser is far 
blue-detuned from a ground state transition, an atom approaching the evanescent field experi- 
ences a repulsive optical dipole force proportional to the intensity of the field: 

C/,«M^,-2-.^ (11) 

4A 

where €Iq is the Rabi oscillation frequency the atom would undergo at the dielectric surface if 
the laser was on resonance, A is the detuning of the field from the transition frequency, and Ky 
is the field amplitude decay rate (typically of order the angular wave number of the laser). For 
sufficiently large Ud, this optical potential may repel incident, cold thermal atoms (e.g. falling in 
vacuum under the influence of gravity) with minimal heating from spontaneous emission. The 
inclusion of van der Waals effects adds an attractive, ~ z^^, potential associated with the red- 
shifting of the atomic ground state (as measured in ll22l ). but also modifies the laser detuning 
as both ground and excited states are perturbed. These surface-induced modifications to the 
optical potential may be used to produce van der Waals-en/ianceii optical potentials, where the 
repulsive energy is derived from the laser fields, but the spatial variation of the potential carries 
the sharp character of the surface interactions: as the van der Waals, ~ z^^, modifications occur 
predominantly in the near field where the un-perturbed Ud is roughly linear, the enhancement 
will be sharply peaked. As a simple example, if the excited state experiences a very large, 
positive van der Waals shift from coupling to resonant surface excitations, the field detuning will 
decrease faster than the red-shifting of the ground state. In the perturbative regime, A ^ 5Ei, 
the van der Waals contributions to the potential are 

f/w/w^^5£; + 5£/, (12) 

and so an appropriately large and positive SE^ will provide a steep (^^ z^^) and repulsive C/,,c/w- 
Also note that the probe need not be evanescent in this case to achieve near field repulsion. 
However, it is important to recognize that this system will also experience an enhanced rate 

of spontaneous emission in the near field w j^^e, which will be comparable to the opti- 
cal potential enhancement if 5Ee ~ TiRe, as suggested in the end of section ID Although this 



necessarily increases the heating rate, surface-induced spontaneous emission enhancement can 
sharply redistribute atomic populations in the near field ifTSl to advantageous effect. For exam- 
ple, enhanced emission may "turn off" target optical dipole transitions and thus quench dipole 
forces. If a ground state atom experiences balanced attractive and repulsive optical dipole forces 
at moderate atom-surface separations 1121 . greatly enhanced spontaneous emission on the "at- 
tractive" dipole transition can produce a net repulsive force in the surface near field, producing a 
sharp barrier and preventing collisions with the surface. As this general approach is particularly 
valid when spontaneous emission enhancement dominants any energy level shifts, its applica- 
tion may be particularly useful when the stability of atomic transitions is required. Moreover, in 
comparison, to other "repulsive" van der Waals interactions that rely on the reflection of atomic 
wavepackacts from steeply attractive van der Waals potentials 1231 . both approaches outlined 
above are designed to manipulate the classical center of mass, and thus may be performed with 
atoms with relatively small de Broglie wavelengths. 

More accurate and general optical forces for thermal atoms may be predicted from the steady 
state solution to the atomic master equation {h— 1) |241 : 



where p is the atomic state density matrix, a„„ is the atomic lowering operator |n)(a|, n„„ 
is the Rabi frequency induced between states a and n, is the energy of the laser photons 
near-resonant with the transition of a ground state atom to level a minus the energy of state a, 
R\"g is the total decay rate of state a to n, and ({•,•}) [•, •] is the (anti-)commutator In the cases 
considered below, the nullspace of the superoperator ^ contains a unique atomic steady state 
Pss, to which an arbitrary initial atomic state will decay in time of order max{/?,7„'}. Typically, 
the internal state of the atom relaxes much faster than the external parameters change (i.e. the 
velocity of the atom times the spatial derivative of the coefficients in Eq. (fTSl l is much less 
than R,ui)- In such regimes, use of Pss to estimate the force on the atom as a function position 
may be justified. To associate an internal atomic state to a force on the center of mass (CM) 
motion of the atom, we adopt a semi-classical approximation that assumes quantized internal 
dynamics (Eq. (fT3]l), but classical CM dynamics. In particular, as the expected force on the 
center of mass of an atom by Ehrenfest's theorem is — Tr([Vi{^j^ ,//]p) (where is the center 
of mass position operator), declaring the center of mass a "classical" variable, r^M^ allows the 
construction of an effective potential for the CM motional dynamics 



This description should be most accurate at capturing the essential contributions to the mean 
force on slowly moving atoms with small de Broglie wavelengths from reactive ("dipole") 
interactions with the laser fields |24il . Similarly, it is straightforward to show that the steady 
state heating rate of the atom due to isotropic spontaneous emission is 



p=if(p) 
Jf(p) = /[p,//]+£^„„(p) 




(13) 




(14) 



3c^kBm 
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where m is the mass of the atom. 
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Fig. 2. Left: diagram of a typical optical atom mirror, as described in section[3] Right: the 
^^K fine structure energy levels relevant to the example in section|4] 

4. Example 

4.1. ^'^KandlTO 

As described in section [2] and appendix |A] resonant van der Waals effects are a product of 
matched, high quality surface polariton and strong optical dipole transition resonances. As is 
suggested by Eq. |20l dipole transition moments tend to increase rapidly with wavelength 
|w„a|^"'). Similarly, compared to plasmonic metals like Au and Ag, relatively stiff materials 
like sapphire, Y2O3, and SiC tend to feature high quality surface phonon-polariton excitations 
near crystal lattice vibrational modes in the far infra-red (<1000cm^') 1251 l26l l27l . Thus, 
strong resonant interactions with common bulk materials are likely to occur in multiply-excited 
states of alkali atoms, which have strong, ^ 10/im wavelength transitions 128] [29l. However, 
working with these highly excited atomic states is often very difficult. For example, although 
the resonant van der Waals interaction of the Cs 6D3 /2 — »■ 1P\ /2 ( 1 2-/im) transition and oriented 
crystal sapphire is very large [131, an even partially excited 6D3/2 state will populate so many 
lower-lying levels through spontaneous emission l28l that construction of an enhanced optical 
potential is difficult both to model and achieve. 

However, there are combinations of common materials and convenient alkali states that admit 
a significant interaction without large complications from the multi-level atomic structure. The 
relatively strong, 3D5/2 ^P^/i 1.178jLim transition in neutral ^'^K 1301 is slightly to the blue of 
a surface polariton resonance recently identified in indium tin oxide (ITO) |[3n[32l . which sug- 
gests that the 3D5 ^2 state may experience a net positive van der Waals energy shift in the vicinity 
of ITO. Due to dipole transition selection rules, 3D^ j2 states spontaneously decay only to 4P3 /2, 
and 4P3 /2 decays only to the /2 ground states [[28 ] . Thus optically probing the 4P3 /2 ^ 451 /2 
and 3D5/2 ^ 4P3/2 transitions forms a closed system between these fine-structure states. More- 
over, cycling transitions between the |45'i/2,^' = 2,m/r = ±2) ^ \\PT^j2^F = 3,m/r = ±3) ^ 
|3D5/2,f = 4,mf = ±4) hyperfine magnetic sublevels are possible (^^K has a 3/2 nuclear spin 
II33I ). Calculation of the van der Waals effects on these fine structure states by a near field pla- 
nar ITO film is done in appendixlAland presented in table|2] The results anticipate a net positive 
energy shift of the 3D5 /2 state, as well as a significantly enhanced spontaneous emission rate 
on the 1.178/im transition. 

In figure [3] we plot the steady state energy expectation (Eq. ( fT4] i) of a van der Waals- 
enhanced optical atom mirror for ^^K above an ITO interface as a function of atom-surface 
separation. Employing a probe configuration that exploits both the van der Waals energy shifts 
and enhanced emission to produce a sharply peaked potential in the surface near field, the 
4P3/2 ^ 451/2 766. 7nm transition is driven with maximum Rabi frequency, detuning, and field 
amplitude decay {ni,Ai, K-i}/27r = {100MHz,50MHz, (767nm)-'} and the 3D5/2 ^ 4P3/2 
1.178/im transition is driven with {i22,A2}/27r = {100MHz,0MHz}. As ITO is quite lossy at 




Fig. 3. Left: effective, steady state optical potentials as a function of K-surface sepa- 
ration for both ITO and (dispersionless) ITO* materials with {fli,Ai, Ki,0.2,A2}/2n = 
{lOOMHz.SOMHz, (767nm)"' , 100MHz,0MHz} and atomic parameters from tables[T]and 
121 Right: steady state atomic level populations for the same. 

1.178/xm (see Eq.l2Tli. the second transition is assumed to be propagating in vacuum (k^ = 0) 
rather than evanescent. The free-space, spontaneous decay rates of these levels are listed in ta- 
ble [T] and the additive enhancement of spontaneous emission rates and energy level shifts from 
van der Waals interactions are listed in table |2] To highlight the significance of the ITO polari- 
ton resonance, figure [3] depicts the effective, steady state optical potential for a ^'^K atom both 
in the near field of an ITO surface and in the near field of a similar, but dispersionless dielectric 
"ITO*" (see appendix lAli. The ITO potential exhibits a greatly enhanced potential barrier over 
the dispersionless surface, capable of deflecting cold atoms falling from mesoscopic heights. 

We begin to understand the origins of the enhanced potential by diagonalizing the "force 
operator" F = for the two potentials as a function of position and considering its eigen- 
values. As plotted in figure H) in its eigenbasis, the diagonal elements of F from the far field 
up to about .Ifim suggest that the two systems feature identical repulsive and attractive op- 
tical dipole force eigenstates, as well as a state that experiences no force. In the near field, 
each force eigenstate attains a van der Waals-like character, with the most significant differ- 
ence between the ITO and ITO* systems emerging from the initially force-less state, due to the 
difference in AEj,o^p calculated for the two surfaces (table |2l): the polariton resonance in ITO 
produces a positive A£'3Dj^, and thus a repulsive F22, while F22 for the dispersionless ITO* has 
a familiar, attractive, van der Waals character. The total expected, steady state forces, and thus 
the effective potentials depicted in figure[3] may be constructed from these position-dependent 
force eigenvalues and the diagonal elements of the steady state atomic density matrix in the 
force eigenbases, Tj,-, . The right half of figure |4] depicts the steady state 77 diagonal elements. 
As to be expected, both the ITO effective potential in figure [3] and the population in the van der 
Waals-repulsed state 7722 in figure |4] begin to decrease with decreasing separation at about the 
same point (~ 80nm). It may appear from figure|4]that the enhancement of the optical forces in 
the ITO system is due to the near-field repulsive F22 eigenvalue. While this is largely true, the 
close correspondence of T]n in the two material systems is a coincidence caused by the large 




Fig. 4. Left: eigenvalues of the force operator as a function of separation for the ITO and 
ITO* systems depicted in firgure[3] Right: steady state atomic populations in the (position- 
dependent) force eigenbasis for the same. 

enhanced decay of the 3D^ j2 state in the ITO system. If this decay enhancement is artificially 
removed in the ITO system, ITO's rjn is significantly smaller over the 60nm-250nm separation 
range, where F\ \ is still repulsive. In this respect, the optical forces in the ITO near field can be 
considered van der Waals-enhanced by both zi£'3o and ^4p, ,,^3d,/t ■ 

We conclude this subsection with a mention of the optical powers required to achieve the 
drives in the above example. It is informative to calculate the saturation intensity for an a ^ « 
electric dipole transition ll34l 

_ hRil(ol2Ja-\-l 

127rc2 2/„+r ^^^^ 

Using the data in tablelT] the saturation intensities of the transitions we are driving are Isai.s^p = 
3.4mWcm^2 and Isai.p^D ~ .47mWcm^2. Using Igat, optical intensity may be related to Rabi 
frequency by / = 2hat{^na / Rm,f ■ Thus, = 27r x lOOMHz requires 1.9Wcm ^ of optical 
power at 767nm and i22 = 27r x lOOMHz requires 618mWcm^2 at 1.178/zm. Although these 
intensities are non-trivial, restricting the illumination region to Imm^^ should bring the power 
requirements comfortably within the range of commercial diode lasers. 

4.2. Quantum trajectory simulations 

After falling 1mm, an atom initially at rest accelerates to .14m/s (for ^^K, this corresponds 
to a kinetic energy of .96MHz). Because the internal state of the ^^K atom should relax at a 
rate of order 24.5/Xs^' (see section [3] and table[T]i and the near field enhancement feature de- 
picted in figure[3]is much broader than the ^lOnm the atom would be expected to travel before 
reaching steady state, we can partially justify the use of the steady state approximation. Sim- 
ilarly, the thermal de Broglie wavelength for a 1/xK (polarization gradient cooled ||35]| ) ^^K 
atomic ensemble is only h{2nmkBT)^^^^ = 7nm, which supports the classical CM approxima- 
tion for a typical cold atom ensemble. On the other hand, as the laser fields in this example 
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Fig. 5. Quantum trajectory simulations of K incident on the optical potentials represented 
in figure[3]after falling from a height of 1mm. The left (right) graph depicts 100 trajectories 
of the atoms falling towards ITO (ITO*) surfaces. The highlighted green and red trajecto- 
ries have the median reflected and absorbed escape times, respectively. The ITO potential 
reflected 88% of the atoms, while the ITO* potential reflected only 7%. 

excite the atoms significantly (see figure O, incident atoms will heat from spontaneous emis- 
sion scattering, complicating use of the steady state approximation. To demonstrate the effects 
of spontaneous emission, illustrate some dynamics connecting the internal and external atomic 
variables, and numerically evaluate the steady state approximation introduced above, we con- 
sider quantum trajectory simulations I36l of ^^K atoms falling from a height of 1mm onto the 
one-dimensional ITO and ITO* optical potentials represented by figure |3] 

Rather than describing the state evolution of a single atom, it is more accurate to interpret 
the atomic master equation ( fTsT i as evolving the mean internal state of an ensemble of inde- 
pendent atoms, without any measurements being made on the system 1371 . This "ensemble 
average" evolution is ill-suited to demonstrate the dynamics between an individual atomic state 
and the forces on the atomic mass. In contrast to its associated master equation, an ensemble of 
quantum trajectories presume that at least partial measurements are being made on the system, 
which add a stochastic element to the dynamics, but helps preserve the state purity of each tra- 
jectory ll36l . In this subsection, we construct such trajectories by simulating experiments with 
photon counting on all radiation fields with perfect efficiency, spectroscopically filtered with 
a bandwidth small enough to perfectly distinguish spontaneous emission from the two atomic 
excited states, but much larger than any of the other dynamical rates in Eq. ( fT3] l. Although 
not experimentally realistic, these trajectory simulations in particular help elucidate the effects 
of spontaneous emission on the motional dynamics of the atoms. As purity is preserved in 
these simulations, the (unnormalized) evolution of a pure, internal atomic state is guided by the 
stochastic Schrodinger equation (/z = 1) ll38l 
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Fig. 6. A .5/xs slice of a representative ITO trajectory. Left: the instantaneous force eigen- 
state populations rj,-,- as a function of time. Right: the expected force on the atom as a 
function of time. Purple and cyan lines indicate different types of photon counting events 
in the trajectory. 

where 1 1/) is the internal state vector of the atom and dQ,ui = lin each time step with probabihty 
^'na {w\ ^na ^na \w)df ™d is zero Otherwise. And as above, we use a semiclassical approximation 
that the evolution of the (classical) atomic center of mass is driven by the force {\if\F {rcM) \ W)^ 
that spontaneous emission events (dQ„a = 1) impart a random momentum kick with a proba- 
bility distribution associated with the photon scattering randomly into An steradians, and that 
rcM determines the dynamical rates in Eq. (T7\ according to the atom-surface separation. 

The K-ITO and -ITO* separation as a function of time in 100 quantum trajectories each 
is plotted in figure |5] Despite an average 114 spontaneous emission events per drop, the ITO 
potential reflects 88% of incident atoms with a turning point that suggests the steady state 
effective potential of figure|3] Similarly, the ITO* potential reflects a minority (7%) of incident 
atoms, as to be expected from the initial kinetic energy (.96MHz) and figure |3] The imperfect 
reflection or transmission of atoms in the two cases is due to "heating" from non-deterministic 
scattering events. 

Figure |6]reveals some of the essential dynamics driving the trajectories, relating the internal 
atomic dynamics and the forces on the atomic center of mass. A .5/is time slice of a repre- 
sentative ITO trajectory is given in both panels (during which the atom-surface separation was 
~200nm). As a function of time, the (instantaneous position-dependent) force eigenstate pop- 
ulations rjii are presented in the left figure. Purple (cyan) lines represent 1.178/im (767nm) 
spontaneous emission events, which instantly project the ^^K atom into the 4P3/2 (451/2) state. 
The right figure depicts the force expectations {y/\F\\if) over the same time interval. Because 
our particular choice of measurement projects the system into an atomic eigenstate, rather than 
an energy eigenstate, the internal states T],-, undergo rapid and (roughly) coherent evolution 
between spontaneous emission events, giving rise to oscillations in the expected force on the 
atoms. Although the exact evolution is specific to the system, separation, and choice of measure- 
ments in the simulation, it is important to note how each photon count abruptly ends the pre- 
vious form of Schrodinger-like evolution and initiates a new form, characteristic of the type of 
count that occurred. Thus we see that in addition to the random momentum kick from photon 
scattering, spontaneous emission also gives rise to fluctuations in the dipole force that, although 
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Fig. 7. Thick lines: quasi-mean r]„ from a single quantum trajectory as a function of sepa- 
ration. Solid (dashed) lines indicate the incident (reflected) portions of the trajectory. Thin 
lines: steady-state rja from figure|4] 

averaged over in the steady state potential Eq. ( fT4b . should be considered another source of heat- 
ing. However, this insight also provides a clear picture of how spontaneous emission rates can 
effect the time-averaged forces on the atom: the total, time-averaged force may be estimated 
from the relative frequency of the spontaneous emission events and the time-averaged force of 
the Schrodinger-like evolution they initiate (over a characteristic time interval). 

Considering whole trajectories again, we re-examine the use of steady state approximations. 
Thick lines in figure |7] depict time-filtered rja as a function of separation in a representative 
trajectory. Low-pass, post-filtered with a cutoff frequency roughly 1/10* the free space spon- 
taneous emission rates, these tj,-, represent the quasi-mean atomic populations relevant to the 
transient, atom mirror interaction. While both incident and reflecting, the quasi-mean tj,-, fol- 
low the underlayed, steady state rj,,- from figure|4](thin lines) with ^^.1 precision. While figure 
|5]suggested the validity of Ueff at a single position (the classical turning point) over many tra- 
jectories, figure|2]suggests the steady state potential for a single trajectory over a wide range of 
separations. 

Finally, we note that other useful, van der Waals-enhanced optical potentials may also be 
constructed. For example, for the same optical parameters as in figure [3] but with an addi- 
tional, -3MHz detuning of the 30^ n ^ ^P^/i transition, a signifiant Ugff minimum forms in 
the near field, as depicted in figure [8] Quantum trajectory simulations of atoms initially at rest 
at the potential minimum suggest a mean trapping time of 20/Xs before heating effects eject 
the atom (figure whereas ITO* potentials produce no such trapping and lose atoms within 
/xs'. Although these trapping times are modest in this (non-optimized) potential, they are long 
compared to many fast atomic processes that characterize near-field atomic systems (e.g. ||39l ). 
Traps and barriers with significantly lower heating rates would be possible for systems with 
anomalous van der Waals interactions that are an order of magnitude or more larger than the 
ground state perturbations (e.g. effects on the order of demonstrated in Cs and sapphire inter- 
actions lITsI ). 

5. Conclusion 

Although encouraging that significant surface-repulsion enhancement appears experimentally 
plausible with anomalous van der Waals effects of only moderate strength (predicted using 
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Fig. 8. Left: effective, steady state optical potentials witli {Q.i.Ai.Ki,Q.2.A2}/27!: = 
{100MHz,50MHz, (767nm)"' , lOOMHz, -3MHz} tliat exhibit a near field potential min- 
imum. Right: 100 quantum trajectory simulations for "''^K atoms initially at rest in the ■'^K- 
ITO potential minimum. The highlighted green and red trajectories have the median re- 
flected and absorbed escape times, respectively. 

published atomic and material data), many of the specific results of the previous section are 
of course limited to the case of ^^K and carefully prepared ITO surfaces f3T\ (although, in 
addition to the surface resonance emphasized here, the unique optical and electrical properties 
of ITO suggest a potential utility in integrated atomic systems in general). Although attractive 
pairings of alkali atoms and common materials are limited by an apparent dearth of high qual- 
ity surface polariton resonances in the near infrared ll26ll . engineering useful resonances with 
even simple surface patterning is a distinct possibility. Guided in particular by Eq. (O, it is 
straightforward to recognize and analyze how construction of a transfer region (figure [1]) with 
planar thin films 1401 , gratings BTl 1421 , or even material dopants 1431 could be used to per- 
turb the dispersion of surface excitations, tuning them to useful atomic transitions. This line of 
research would naturally compliment current investigations that leverage similar effects using 
plasmonics to enhance and collect the spontaneous emission from embedded dipoles Il45ll44l . 

Harnessing surface effects for strong, near field manipulation of gas-phase atoms would en- 
able integrated atomic systems and (at least in this implementation) employ an interesting in- 
tersection of fundamental QED, atomic physics, material science and surface polaritonics. Al- 
though exceptional, anomalous van der Waals effects admit an intuitive physical picture (sec- 
tion |2|i that should facilitate their generalization. And as suggested by section H] the effects 
need not be large for significant enhancement of near field potentials, and efficient steady state 
analysis may be reasonably accurate. However, we emphasize again that the identification of a 
broad class of surfaces with resonant excitations tunable to atomic transitions would open the 
door to general application of van der Waals-enhanced optical dipole potentials for atoms. 

Appendix 

A. Calculation of van der Waals effects 

Estimating the influence of a surface on atomic states is complicated, even in the near field 
where the field susceptibilities may be approximated by Eq. In this limit, with a simple 



planar interface, the energy shifts and spontaneous emission enhancement are 
oEa = -> -— ^ -/ dQ +2Re— — — 0(a)„„) 

^ -^m„„(a;;{+a,';jz-3 ds) 

n 

^ + iMff (^»») - 1 Q( ) ^ rz-^M„„r„„z-^ (19) 

16/!z3 ei((o„„) + l 

We can therefore break up our calculation into separate parts: the dipole transition strengths, 
Man, and the image factors, t^na and Vna- 

We will neglect any atomic hyperfine structure and assume an isotropically prepared atomic 
dipole for calculation and conceptual clarity. The widths of e(a)) features will be on the order of 
.leV, while typical alkali hyperfine splittings are only of order lOGHz or less 1331 . Thus image 
factors should not vary significantly between hyperfine states. However, these splittings are 
of order the energy scales of typical laser Rabi frequencies and detunings. Moreover, atomic 
dipole moments vary greatly between different hyperfine spin projections. A more thorough 
analysis of Ma„ should be considered in any experimental implementation. 

Under these approximations, we may estimate the atomic dipole strength with 

,™|2 , |,,fl«|2 _ ftp/i- ( C '- -Jn-Jn 



+ (l+2^0(co„.)) (20) 

where /jj^ is the free space decay rate for either a ^ « or a <— n, and 7,- is the spin-orbit 
quantum number of state /. Using atomic data from 1301 and l28l . we list the calculated Mq„ for 
the 45i/2, 4^3/2, and 3D5/2 levels of neutral '^'^K in table[l] 

The image factors may likewise be calculated from the bulk material's dielectric function. In 
modeling the response of indium tin oxide (ITO), we follow 1311 and use an empirical Drude 
free-electron model 

e,{(o) = E^- ' . . (21) 
ct)(a) + i7) 

with measured material constants = 3.8, (Op = 2.19eV, and J— .llleV. Using Eq. (1211) . we 
list the calculated A„„ and r„„ in table [T] 



Table 1. Calculation of the transition dipole strengths and image factors for the 45i/2i ^^3/2' 
and 31)5/2 levels of neutral ^'^K and ITO. Data taken from ll30ll28il3Tl . 



a n 




Man (kHzjUm^) 




45i/2 4Pi/2 

4^3/2 
5A/2 
5i^3/2 


1.6093 36.9 
1.6165 37.2 
3.0613 1.98 
3.0637 2.14 


.6758 
1.3447 
.0053 
.0114 


.7373 

.7368 
.6821 
.6821 


4P3/2 45i/2 

5^1/2 

303/2 
305/2 


-1.6165 37.2 
.9894 14.2 
1.0527 4.09 
1.0524 24.5 


.6723 
.5595 
.2627 
2.4067 


-.7368 .6570 .1133 

.7856 

.7793 

.7794 


3D5/2 4P3/2 

5^3/2 


-1.0524 24.5 
.3948 1.58 


1.604 
1.3068 


-.7794 -2.3873 7.6597 
.8741 



Table 2. Total van der Waals energy shifts and enhanced decay rates for K, ITO, and 
ITO* as determined by Eqs. J18t and U9b . and table[T] 

dEa,Rm, ITO (kHzjUm^) ITO* (kHzjUm^) 





-1.5 


-1.2 




-2.3 


-2.3 




3.9 


-1.7 




27rx .76 







27rx 12.3 






The aggregate van der Waals effect for each atomic level is determined by Eqs. (fTSI l and ( fT9] l, 
and the entries in table [T] The results are presented in table |2]in comparison to a dispersionless 
dielectric called "ITO*" with ei (o) = £„. The largest effect is the enhanced decay R4p^^^^j,o^^^ . 
The nearby surface polariton resonance in ITO (ei (1.0524eV) = —0.4826 + 0.4517;) and the 
relatively strong 3D5/2 ^ '^P^/i dipole moment (see table[T]i combine for a net positive energy 
shift and a significantly enhanced spontaneous emission rate. 
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